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Рассмотрены вопросы выбора и построения вычислительных ресурсов и способы эффективного их 
использования при приближенном вычислении интегралов от быстроосциллирующих функций 
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Цель работы. Поскольку «легкие» задачи уже перерешали и остались «трудные» 
задачи, которые невозможно решить при помощи штатного материального обеспече- 
ния с требуемым качеством, необходимо привлекать для ее решения резервы оптимиза- 
ции вычислений [1]. Как их использовать и в каком порядке, подсказывает компьютерная 
технология решения задач с заданными значениями характеристик качества. 
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Постановка задачи. Ниже будут приведены элементы компьютерной технологии 
построения вычислительных алгоритмов (в.а.) с заданными значениями характеристик 
качества для решения задачи интегрирования быстроосциллирующих функций вида: 

—@х 
Ь 
Ко) = | 1(х)4 5т ах ах, 


ы с05 @х 


(1) 
и =5 | а 
= И Хх, 
где /(х)ЕЁ (РЁ — некоторый класс функций); @ — произвольное вещественное 
число (@ >2л(Ь-а)). Информация о /(х) задана № значениями в узловых точках ф | в. 


из области ее определения: р = Г; т =; характеризует точность задания /(х, ): 


Я-л 


<=; 1=0,№М-1; и(х) — вейвлет-фукнция [2], которая колеблется около оси х, удов- 


[.) 
летворяет условию 7 х ах = 0 и быстро стремится к нулю на бесконечности. 
—с 
Пусть с(У) — модель компьютера, с(У)еЕС(У) (С(У) - класс моделей компью- 
теров). 
Общая ситуация построения вычисления интегралов (1) с точностью &,(#>0) 
при ограниченных вычислительных ресурсах описывается условиями 


ЕОХ.У) <. (2) 
и (3) 
М(Г,Х,У)<Мо(=), (4) 


где =,Ту,Мо — заданные числа; Е(1,Х,У) — полная погрешность приближенного 
вычисления интегралов (1), которая является суммой трех составляющих: Ё, (+) — неустра- 
нимой погрешности за счет неточности входных данных; Ем(-) — погрешности метода; 
Е+(.) — погрешности округлений; Х,У -— векторы параметров, которые характеризуют 
соответственно в.а. и компьютеры из классов А и С; Т(Е,1,Х,У) и М(Е=,1,Х,У) - 
процессорное время и память компьютера, необходимые для приближенного вычи- 
сления интеграла. 

Приближенное значение интеграла (1), для которого выполняется условие (2), 
называется 5 -решением, А(=,Х,У) — множество в.а. построения = -решения в данной 
модели вычислений. 

Вычислительный алгоритм, который удовлетворяет условиям (2), (3), будем назы- 
вать Т-эффективным. 

Далее будем предполагать (если не допускается другое), что память М может быть 
расширена до необходимого объема, то есть ограничение (4) можно снять, но, скорее 
всего, за счет увеличения характеристики Т (процессорного времени). Это может быть 
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осуществлено, например, путем наращивания части «медленной» дисковой памяти 
в общей структуре компьютерной памяти. Учитывая тот факт, что при =—>0, Му о 


(например, когда идет речь о погрешности округлений или погрешности метода), будем 
считать, что #250 >0, где =0 — задано. 


На практике используют некоторые оценки полной погрешности в.а., ее составляю- 
щих и процессорного времени. При этом различают оценки априорные и апостериорные, 
мажорантные и асимптотические, детерминированные и вероятностные. Возможность 
и целесообразность использования указанных оценок и способов их получения зави- 
сит от типа, структуры и точности априорных данных задачи и в.а. от того, с какой целью 
вычисляется оценка, а также от вычислительных ресурсов. 

Мажорантные априорные оценки гарантируют верхнюю грань оцениваемой вели- 
чины и выражаются через известные величины, вычисление их не требуют значительных 
вычислительных затрат, но значения оценок часто очень завышены, поэтому выводы 
на их основании касательно возможности вычисления интеграла при условиях (2), (3) мо- 
гуг быть ошибочными. 

Асимптотические оценки аппроксимируют величину, которая оценивается. Варьи- 
рованием параметра может быть достигнута желаемая близость оценки к оцениваемой 
величине, но вычисление таких оценок связано со значительными вычислительными 
затратами. Эти оценки, как правило, апостериорные. 

В алгоритмическом обеспечении вычисления интеграла (1) при соблюдении усло- 
вий (2), (3), в соответствии с указанными характеристиками оценок, нужно предусмо- 
треть возможности вычисления разных видов оценок характеристик Ё(ЁЕ„,Ем,Е.) [3]. 


При ослабленных ограничениях (2), (3) достаточными оценками могут быть менее точные 
и менее сложные (в вычислительном плане). При усиленных ограничениях (2), (3) 
используются асимптотические (апостериорные) оценки. Например, условие (3) может 
привести к высоким требованиям к точности параметров вычислительного процесса, 
которые вычисляются на основе оценок погрешности решения. 

Квадратурные формулы вычисления интегралов типа (1) широко используются 
при решении следующих важных классов прикладных задач: прикладной статистики, 
анализа сейсмограмм землетрясений, анализа и синтеза речевых сигналов, решения 
краевых задач для уравнений в частных производных, информационной безопасно- 
сти и др. Поэтому приближенное вычисление интегралов (1) для Г(х) из некоторых 
классов и анализ точности соответствующих квадратурных и кубатурных формул 
являются актуальными задачами вычислительной и прикладной математики. 


Компьютерная технология 


В случае, если не удается построить в.а. ае А(=,Ту) вычисления = -решения 
задачи (1) на классах функций ГеР, Ру, Еу; [4], для окончательного вывода о 


возможности такого решения важно иметь точные оценки снизу (или близкие к ним) 
приближенного решения и вычислительной сложности задачи. Воспользовавшись 
этими оценками в.а., можно сделать окончательный вывод: решение задачи с задан- 
ными значениями характеристик качества можно построить, или что такое решение 
построить невозможно и, возможно, необходимо /(х) погрузить в более узкий класс 


функций (например, Ру или Ру -) или использовать компьютер другого класса. 
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Для улучшения качества квадратурных формул вычисления интегралов от быстро- 
осциллирующих функций (1) целесообразно использовать следующие резервы опти- 
мизации вычислений: 

1. Резервы уменьшения меры погрешности: 

а) за счет уточнения входных данных (ЕЁ, ); 


— уточнение класса подынтегральных функций; 

— корректировка входной информации; 

— повышение точности входной информации. 

6) Метода Ем: 

— использование оптимальных и близких к ним по точности в.а.; 

— оптимизация информационного набора функционалов (например, оптимальная 
сетка узлов); 

— увеличение количества функционалов в информационном наборе; 

— использование иного информационного оператора; 

— полное использование входной информации для сужения класса подынтеграль- 
ных функций. 

в) Округлений (Е): 

— использование методов оптимального порядка точности для высокоточных вычи- 
слений; 

— сужение классов (с целью уменьшения количества функционалов в информа- 
ционном наборе); 

— использование схем вычислений, которые минимизируют скорость накоп- 
ления погрешности округлений; 

— увеличение оптимального порядка точности путем выбора класса информацион- 
ных операторов; 

— выбор и моделирование правила округления. 

2. Резервы уменьшения процессорного времени: 

— увеличение длины разрядной сетки; 

— использование оптимальных и оптимальных по порядку точности квадратур- 
ных формул; 

— использование уточненной входной информации (как следствие — менее жесткие 
ограничения на погрешности метода и округлений); 

— улучшение качества оценок погрешностей метода и округлений; 

— улучшение точности вычисления параметров вычислительного процесса; 

— использование алгоритмов быстрых ортогональных преобразований; 

— использование «быстрой» арифметики; 

— согласование в.а. с архитектурой компьютера; 

— распараллеливание вычислений. 

Необходимость в вычислении интегралов (1) с необходимыми значениями характе- 
ристик качества и в получении информации о возможности обеспечения такого реше- 
ния возникает, например, в следующих случаях: 

— необходима диагностика качества полученного значения интеграла по точности, 
времени и необходимой памяти компьютера; 

— необходимо до нахождения значения интеграла знать о возможности его нахож- 
дения с необходимыми ограничениями на Ё, Т, М ; 


— необходимо обеспечить необходимую точность вычисления интеграла (1). 
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Программа для вычисления значения интеграла может принадлежать к одному 
из следующих типов: 

1) программа нахождения приближенного значения интеграла с заданной точ- 
ностью. Значения ее управляющих параметров, которые соответствуют необходимой 
точности = > 0, вычисляются самой программой; 

2) программа, которая вычисляет = -решение задачи и вычисляет оценки погреш- 
ности Е(1,Х,У) полученного решения; 


3) программа, которая находит приближенное (с точностью &> 0) значение 
интеграла (без вычисления оценки погрешности). 
Каждая из приведенных программ может вычислять апостериорные оценки Е Т; М. 


Для того, чтобы еще до решения задачи можно было судить о возможном обеспе- 
чении ее решения с необходимыми характеристиками качества, необходимо иметь соот- 
ветствующую программу вычисления априорных оценок ее характеристик Ё, Г, М, 
оценок наилучших управляющих параметров Х’ (если такие у нее есть), а также выдачи 
информации о значениях оценок ее характеристик, полученных в результате тестирования. 

Технология вычисления интеграла обеспечивается пользователем в диалоге 
с компьютером с помощью программ решения задачи и вычисления априорных оценок 
ее характеристик. Она существенно зависит от типа программ, к которому она принад- 
лежит, а также от наличия программы вычисления оценок характеристик. 

Так, если программа вычисления интеграла принадлежит к первому типу, то 
при заданной (из допустимого диапазона) точности вычисления интеграла вычисляют- 
ся оценки ТГ и М . Если они соответствуют заданным ограничениям (3), (4), интеграл 
может быть вычислен с заданными значениями характеристик качества. 

В случае, когда программа принадлежит ко второму или третьему типу, то на осно- 
вании вычисленных априорных оценок А, Т, М можно судить о возможности выполнения 


ограничений (2) - (4). Если программа не вычисляет априорную оценку Ё ‚ а только 
оценки 7 и М, то вычисляются только они. Если ограничения на Т и М выполняются, 
то задачу можно решать выбранной программой. При этом для программы второго типа 
при помощи вычисленной апостериорной оценки точности можно сделать вывод о вы- 
полнении ограничения (2). 

О выполнении (для данной программы) ограничений (2) - (4) можно сделать 
вывод и при помощи результатов ее тестирования [5]. 

Для обеспечения = -решения задачи (1) с заданными значениями характеристик 
качества (2) — (4) необходимо знать, что нужно делать, как это делать и в какой после- 
довательности. Приведем описание последовательности шагов технологической схемы, 
которая реализует этапы приведенной выше схемы решения задач. 

Шаг 1. Выбирается интеграл (см. (1)), который необходимо приближенно вычи- 
слить. Формулируются требования к = -решению. Фиксируется класс подынтеграль- 
ных функций (К, Еу, Ем )- 

В дальнейшем (при возвращении к шагу 1) может быть рассмотрен более узкий 
класс подынтегральных функций с целью расширения возможности удовлетворения 
условиям (2) - (4). 

Шаг 2. Находится (или улучшается) оценка меры погрешности искомого # -реше- 
ния за счет неточности входных данных ( ЕЁ’, ), используя компьютерную модель вычи- 
слений, входные данные и оценки погрешности ее задания. 

Шаг 3. Для найденной оценки Е, проверяется выполнение условия 


Е, <а; =, 0<а, <1 (например, @, =13). (5) 
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Если условие (5) выполняется, перейти на шаг 5, иначе — на шаг 4. 
Шаг 4. Для обеспечения & -решения необходимо изменить (если это возможно) 


требование на точность (увеличить => 0) или улучшить (уменьшить Ё’), используя 
определенные резервы. Перейти на шаг 2. 
В случае, когда Е„ больше заданного @, = и не может быть уменьшено, необ- 


ходимо изменить постановку задачи, то есть перейти на шаг 1, или прекратить поиск 
= -решения в связи с невозможностью его вычисления в данной компьютерной модели 
вычислений. 

Шаг 5. Изучить вычислительную ситуацию для решения задачи. 

5.1 Для вычисления интеграла неизвестен в.а. и соответственно теоретически 
не изучены оценки его характеристик. Перейти на шаг 6. 

5.2 Известен метод вычисления интеграла, но не построены на его основе в.а. 
и программы для соответствующего компьютера с(У) и теоретически не изучены 
оценки их характеристик ( Е, Т, М ); перейти на шаг 7. 

5.3 Для вычисления интеграла имеются алгоритмы — программы из некоторых 
известных программных систем (например, Маетайса, МаНаб, Маса4), но они 
не обеспечены оценками характеристик, за счет использования которых можно было 
бы сделать вывод о возможности найти & -решение интеграла, удовлетворяющее усло- 
виям (2) - (4); перейти на шаг 13. 

5.4 Для вычисления интеграла имеется пакет прикладных программ (см., напри- 
мер, [3], в котором есть в.а. — программа нахождения = -решения задачи; особенность 
этого пакета состоит в том, что его в.а. — программа нахождения = -решения задачи 
сопровождается оценками характеристик Ё, Т, М ; если она обеспечивает нахождение 


& -решения с заданной (в некотором диапазоне О/=) точностью &ЕД\(=), то перейти 
на шаг 14, иначе — на шаг 7. 

Шаг 6. Для нашей задачи не разработан метод ее решения, и его необходимо 
создать. Перейти на шаг 7. 

Шаг 7. Для решения задачи вычисления интеграла известны (или разработаны 
в соответствии с шагом 6 методы (или метод). 

Из множества известных методов приближенного вычисления интегралов, исполь- 
зуя известные априорные оценки их характеристик (например, оценки погрешности 
метода, сложности), выбирается «лучший» по точности и быстродействию, который 
может быть использован для вычисления интеграла. С этой целью используются извест- 
ные знания о классе задач, которые имеются в компьютерных базах знаний, статьях, 
обзорах, монографиях (см., например, [4], [6]). 

Шаг 8. Выбираются параметры компьютера, на котором будет решаться задача, 
исходя из начального анализа ее сложности и требований к характеристикам качества 
= -решения, то есть, доопределяется компьютерная модель с(У)еЕС(У) путем выбора 


параметров У : количество процессоров (А =1 или А>1) и их тип, длина машинного 
слова, количество машинных слов, которые используются при реализации арифме- 
тических операций, оценки времени выполнения операций, режим вычислений, огра- 
ничения Мо на используемую оперативную память. 

Шаг 9. На основе разработанного или выбранного метода разрабатывается 7-эф- 
фективный в.а. (см. [7]) вычисления = -решения задачи и при этом находятся оценки 
его характеристик при помощи следующей технологической схемы. 

9.1 На основании выбранного метода для модели компьютера с(У) разрабатывает- 
ся некоторый в.а. (с той или иной детализацией его описания, например, пошаговое 
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описание, описание на алгоритмическом языке), или разрабатывается новый или улуч- 
шается разработанный (при повторных реализациях шага 9.1) вычислительный алгоритм. 

9.2 Находится или улучшается известная априорная оценка М(=,1,Х,У) необходи- 
мой оперативной памяти компьютера, которая будет использоваться в.а. при построении 
приближенного решения задачи. 

9.3 Проверяется выполнение условия (4). Если оно выполняется, то перейти 
на шаг 9.4; иначе — на шаг 8 (или 9, если алгоритм фиксирован). Это делается с целью 
выполнения условия (4) за счет модификации разработанного или построения нового в.а., 
возможно с выбором соответствующего компьютера с(У). 


9.4 Находятся или улучшаются априорные и (или) апостериорные оценки: ЁЕ,, Ех, 


Ех. - соответственно за счет неточности входной информации, метода (алгоритма), 


округлений — при использовании в.а. ае А и модели компьютера с(У). (Примеры по- 
лучения таких оценок в [4], [6]. 

9.5 На основании оценки Е проверяется, выполняется ли условие (2), то есть, 
может ли быть вычислен интеграл с точностью = при заданной информации /,. Если 


= -решение не обеспечивается, то выясняются причины этого и обсуждаются способы 
их устранения (осуществляется анализ, начиная с шага 7 или 9.1 данной технологии). 
При обеспечении нахождения = -решения выбранным в.а. перейти на шаг 9.6, в против- 
ном случае поиск в.а. прекращается и указывается, почему невозможно обеспечить 
выполнение условия (2). 

9.6 Находится или улучшается априорная оценка процессорного времени Т(=,[, ХУ) 
в.а. (Примеры таких оценок см., например, в [7]). 

9.7 Используя оценку Т(=,1,Х,У) в.а., выясняется, обеспечивает ли этот в.а. = -ре- 


шение задачи при условии (3). Если & -решение обеспечивается при этом условии, то тем 
самым будет установлено, что найденный в.а. принадлежит множеству (1) (=,10) е А(Е) 


эффективных по быстродействию (7-эффективных) в.а., которое зависит от выбора Ту. 


Если же в.а. не является Т-эффективным, то выясняются возможности обеспечения 
его Т-эффективности за счет использования резервов оптимизации в.а. (перейти 
на шаг 9.1 или 9.6). 


Замечание 


Если, используя приведенную технологию, не удается построить в.а. вычисления 
= -решения, для окончательного вывода о возможности такого построения важно иметь 
точные оценки снизу (или близкие к ним) погрешности приближенного решения и вы- 
числительной сложности задачи вычисления интеграла. Используя эти оценки в.а., 
можно сделать окончательный вывод: & -решение можно построить или такое решение 
построить невозможно при данных условиях вычислений и, скорее всего, необходимо 
использовать более точную исходную информацию, большего объема или использовать 
другой компьютер. 

Если разработан Т-эффективный в.а. с указанными выше значениями характеристик, 
то перейти на шаг 10, иначе — на шаг 9. 

Шаг 10. Разработанный Т-эффективный в.а. реализуется в в.а.-программе на со- 
ответствующем алгоритмическом языке для выбранной модели компьютера. 

Разработанная в.а.-программа, которая обеспечивает решение задачи с необхо- 
димыми значениями характеристик качества (2) - (4), должна принадлежать к одному 
из трех указанных выше типов программ. 
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Шаг 11. Осуществляется тестирование разработанной в.а.-программы и построен- 
ных для нее оценок характеристик Ё и Г. Тестирование осуществляется на соответству- 
ющем наборе тестовых задач и в соответствии с технологией тестирования характеристик 
качества в.а.-программы (см., например, [5]). 

Если результаты тестирования подтверждают тот факт, что разработанная (или вы- 
бранная) в.а.-программа может обеспечить построение = -решения задачи с необходимы- 
ми значениями характеристик качества (2) - (4), то перейти на шаг 15, иначе — на шаг 12. 

Шаг 12. Проводится анализ в.а.-программы с целью выявления резервов ее улуч- 
шения (оптимизации) по тем характеристикам, которые не удовлетворяют заданным 
условиям; осуществляется оптимизация в.а.-программы и для полученной модификации 
выполняются действия, указанные на 9-м шаге; если исчерпаны все резервы оптимизации 
для выбранного в.а. и соответствующей программы, но ее характеристики не удовлетво- 
ряют заданным ограничениям (2) - (4), то необходимо продолжить поиск необходи- 
мого в.а. в соответствии с предыдущими шагами (то есть, переход на шаг 8 или 9). 

Шаг 13. Для выбранной программы вычисления интеграла проводится (если есть 
необходимость) ее тестирование и получение оценок характеристик: Е — точности 
вычисления интеграла, Г’ времени для вычисления интеграла в соответствии с техно- 
логией тестирования (см., например, [5]). 

Если выбранная в.а.-программа в соответствии с результатами тестирования может 
обеспечить построение = -решения с необходимыми значениями характеристик качества 
по точности и быстродействию, то перейти на шаг 15, иначе — на шаг 5.2 или 5.3. 

Шаг 14. Используя оценки характеристик Е, Т, М в.а.-программы вычисления 


интеграла с точностью = , проверяется, может ли она обеспечить & -решение с задан- 
ной (2) точностью &ЕД(=) и при заданном ограничении на компьютерное время (3). 


Если эта в.а.-программа может обеспечить вычисление & -решения задачи, то перейти 
на шаг 15, иначе — на шаг 5. 

Шаг 15. Вычисление (разработанной или выбранной в.а.-программой) = -реше- 
ния с заданными значениями характеристик качества (2) - (4). 

Шаг 16. Интерпретация результатов вычисления. 


Выводы 


Предложенная и пошагово описанная информационно-компьютерная технология 
вычисления интеграла с заданными значениями характеристик качества решения (2) -— (4) 
дает возможность: 

— осуществить системный подход как к постановке задачи, так и к ее решению; 

— изучить вычислительную ситуацию (см. шаг 5) при вычислении интеграла и 
выявить все шаги, которые необходимо осуществить, чтобы обеспечить решение за- 
дачи с заданными значениями (2) - (4) характеристик качества решения; 

— при данной вычислительной ситуации осуществить разработку соответствую- 
щей в.а.-программы для решения задачи; 

— использовать соответствующую программу из известных программных систем 
(например, Маетайса, МаНаБ, Ма саа), осуществляя (по необходимости) тестиро- 
вание для выявления той программы, которая обеспечивает решение задачи с задан- 
ными значениями характеристик качества решения; 

— использовать предложенную информационно-компьютерную технологию при ре- 
шении задачи с повышенными требованиями к значениям (2) -— (4) характеристик ка- 
чества ее решения. 
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